Absence of resonant enhancements in some inclusive rates 
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A toy model is defined and solved perturbatively with the aim of examining some claimed "reso- 
nant" enhancements of certain reaction rates that enter popular models of leptogenesis. We find: a) 
that such enhancements are absent; and b) that the perturbative solution, as done correctly using 
finite-temperature field theory is well defined without the "resumming" procedures found in the 
literature. The pathologies that led to the perceived need for these procedures are an artifact of 
uncritical use of weighted vacuum cross-sections in the determination of rates, without adequate 
attention to the effects of the medium upon the single particle states within it. 
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THE MODEL AND ITS SOLUTION 

M 

In some approaches to leptogenesis calculations based on the assumed existence of heavy leptons, certain "resonant" 
effects play a key role |lj-|l(Tj. In the present note we present a toy model in which, we believe, one can examine 
authoritatively some assumptions that are often made in the dynamical treatment of these models, but which we shall 
find are false. We also intend our results to serve as a cautionary warning against the uncritical calculation of reaction 
rates in media by averaging vacuum cross-sections, calculated perturbatively, over unperturbed thermal distributions 
00 . for the participating particles. 
V} ■ We consider a set of particles described by: 



CN \ 1. A scalar (Higgs) field, 

■ 2. A heavy neutral lepton field N. 

3. A light neutrino field L. 

4. A light "anti-neutrino" field L. 

We can make our dynamical point most economically in a non-relativistic framework in which the fields N,L,L all 
i-C ■ contain only annihilation (positive frequency) terms. If we were to constructing sums of dangerous Feyman graphs, 
, or those which enter the pathologies that we will be discussing, there would be a detailed correspondence between the 
k>( ' graphs produced in the models of the literature and those produced by the model below, where the particle annihilated 
rS | by L is indeed the anti-particle of the particle annihilated by L. 

The couplings that we shall include are those that enable the processes: N «-> L + $ with a strength parameter g 
that is small enough to allow perturbative development in powers of g; and A" <-> L + $ with a strength parameter g' 
which is taken as very small. We seek results for the rate of change of the number of Z's, to lowest order g' 2 . The 
interesting question is now how the interactions N <-» L + $, with strength g perturb the lowest order result; the 
usual "resonant" calculations find terms that are effectively of lower order in g than the superficial estimate g l2 g 2 . 1 
The physical context will be the following: without tracing the system back to its creation, I take it to have evolved 
long enough, and to be subject to external influences that change slowly enough, for the N, L, $ system to have 
reached equilibrium. In contrast, the L sector is not in equilibrium; for simplicity we take the number of L's to be 
zero at our initial time. The qualitative problem to be addressed originates in the fact that there are reactions like 
L + <f>^L + <I>in which the intermediate N is unstable. This raises some familiar questions about the treatment 
of unstable particles in field theory. But we would prefer to call them "familiar pseudo-questions", in view of our 
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1 To make clear the specific correspondence of this remark with leptogenesis scenarios in the literature, we cite, e.g., two equations from 
Ref. 0; Eq.(13), in which superficial dependence of a relevant averaged cross-section, is proportional to a factor (Ji<h) 2 , where h is a 
coupling constant, and Eq. (18) where the reaction rate (or the leading term thereof) is found to be of order h^h. Equivalent equations 
occur in most of the references that we cite from this subject area. 



2 



contention that thermal field theory turns its crank in the same way, oblivious to whether particles are stable or not, 
as long as we calculate rates that are sufficiently inclusive. 

To exhibit the general structure of such processes requires only one particle in each reaction to have a distribution 
of energy in its states. We take an L or an L of momentum pj to have energy Ej. The L's will taken in a Fermi 
distribution, np(ej) = [1 + exp(/3ej)] _1 where ej = Ej — fi^. We take the $ particles all to have a common energy, 
E$. Taking the number of $ states in a unit volume to be unity, the equilibrium number density of $'s is just the 
Bose function = [exp(/3e$) — I] -1 where e$ = E$ — We take the N particles to have the same number of 

states, but with the Fermi distribution, n^(ejv), where ejy = En — pn- Finally we take the species N, L to be in 
chemical equilibrium, pn = + PL- The L's will remain out of equilibrium. 

For the L and L fields we take the interaction picture representations, 



j 

L(x,t) = {Vol)- 1/2 J2^ e ' iEit ^ Py> 



(1) 



where Cj and Cj are the respective annihilation operators for L and L. Representing N and <I>, with annihilation 
operators, and &i, respectively we write the free Hamiltonian as, 



H = E N J2 4 a i + E *J2 b l b * + E o ( c h + etc 



(2) 



A local interaction of the fields L and L with the fields N and <&, where the latter two describe particles situated 
at the origin, is given by, 



\f(Ej) 



lp(E 3 ) 



1/2 



(c 3 +ec 3 )[N^] 3 +gY, 



lp(E 3 



1/2 



(3) 



where we have defined, 



[N(ty = / d 3 x e- ip J- x A^ t (x,t)$(x,i) (4) 

Jv 

[JV(i)4 f (t)]j = / d 3 x e lp J x A^(x,t)$ t (x,t) (5) 
Jv 

The function f{Ej) in J2Jl includes the density-of-states factor p(Ej), defined such that J2k v(^k) = / dEk p(Ek)y{Ek). 
It includes as well as any momentum dependence of the effective coupling constant, and perhaps an ultra-violet cut-off 
factor. 2 We choose parameters such that En — E<s> > Ej for some range of L, L states. The coupling that is "lepton 
number violating" in the sense that we defined above is given by g' = eg. 

In our detailed example, in order to address a problem that is absolutely well defined, we take the N, <&, L to 
be in thermal and chemical equilibrium. We define pn,P$>,Pl as the chemical potentials for these species, with 
Pn = P<s> + Pl, and we define the quantities, epj — En — Pn and e$ = — In the discussion section we shall 
return to the relationship between the conditions assumed above and the somewhat different conditions that obtain 
in a typical leptogenesis calculation. But we emphasize that our point is not to recapitulate these calculations, but 
rather to address in the simplest terms the dynamical questions raised below. 

Now we easily see the problem of the Feynman graph approach. To first order g' 2 and in the zero'th order in g 2 
the rate of production of L's per unit volume, where the initial state has no L 's is, 

w (0) = 27rg' 2 n F (e N )[l + n B (e*)]f(E N - E 9 ), (6) 



2 Of course all local, relativistic models will have vacuum processes that give ultraviolet divergences that need to be removed by renor- 
malization. But these are irrelevant to the issues we deal with here. 
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where np is the Fermi distribution, np(e) — [1 + exp(/3e)] 1 , ub is the Bose distribution, ns(e) 
and j3 = 1/T. To first order in g 2 we have two processes that produce L's: 



[exp(/3e) - 1] 



A. The process, $ + L^N^$ + L with the rate 



w 



(1) 



2irg 2 g' 2 n B (e.s,)[l + ns(e*)] / dE, 



npjE.-HL) f{Ei)< 
' (E N ~ E* - E t ) 2 



(7) 



B. The process, N + L^Q + L + L^N + L. For pedagogical simplicity we arrange the initial densities such 
that this process is insignificant compared to that of A. We point out in the subsequent treatment exactly the way in 
which the perturbative anomalies in the rate for this second process get cancelled. 

The prevailing view in the literature is that the singularity in the integrand in J2J is regulated by a "resummation" 
that makes the replacement, 



1 



{E N -E*- Ei) 



1 



E N - E 9 - Ei 



where 7 is the decay width of the N particle, followed by another replacement, 



1 



En — E$ — Ei — i-f 



E 



N 



E& — Ei 



{E N -Es,-E t ) 2 + 1 2 



(E 



N 



Eib — Ei 



{E N -E*- E t ) 2 + 7 2 



(8) 



(9) 



The first replacement, above, is the obvious one to take into account the nature of an unstable particle in an inter- 
mediate state. The second replacement is supposed to remove a perceived double-counting, where a term in the the 
resonance production Ti in the reaction L + $ — > L + $ has recapitulated the decay rate from the direct decay of an 
N. We compute 7 from the decay rate for N — > L + $ alone, and it is of order g 2 , whence, 



f , P ^. KEinEN-Ey-Ei) 2 2 

J dEi nF{Ei) [(E N - 23. - Ei f + 7 T = 0(9 ] (10) 

With this modification, of Q is of order g' 2 , or reduced by two orders of g from the superficial behavior of the 
graph, by virtue of the singularity regulation. We do not elaborate further on this procedure, which we claim to be a 
plausible series of ansatzes, but incorrect. 

To find the correct expansion in powers of g of the complete rate calculated to order g' 2 , we show first that the field 
equations themselves give a simple formal expression, valid to order g' 2 , for the net rate of change of the number of 
.L's (the analogue in our model of the rate of the change of lepton number in the theories which we are explicating). 
This expression is of the form of a Fourier transform of an expectation value of a product two "currents" evaluated 
at different space-time points, and is valid for any configuration of the bystanding medium in which we are to take 
the expectation value. 

We define Q = J2i an d n °t e that Q would be conserved but for the eg = g' terms in H. We have 



^=i[H',Q]=i(K 1 -K 2 ), (11) 



where 
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We wish to calculate (dQ/dt) where the bracket indicates the average value in the medium. In the absence of the 
coupling g' we choose the initial medium to be in an eigenstate of Q (we shall ultimately take Q to be zero, as in the 
examples of and J7J), but this is not necessary.) Then the expectation values of K x and K 2 are zero. The terms 
of order g' 2 in the expectation value for dQ/dt are generated by considering the linear response to the part of the 
Hamiltonian that is proportional to g' , namely, SH = K\ + K 2 . 

The linear response of a Heisenberg operator A(t), evaluated at t — 0, to a perturbation in the Hamiltonian of 5H 
is given by, 

5A(Q) = -i [ dt[A(0),SH(t)]. (13) 
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Thus we can exhibit the term that is second order in g' in the expression for the operator dQ/dt, as, 



dQ/dt = 



dt 



(x 1 (0)-/f 2 (0)),(^ 1 (t) + ^ 2 (0 



(14) 



Following the steps of Ref. .11], appendix A, and noting that the terms involving two Ki's or two -ftVs in l)14|l each 
give vanishing contributions, we use time translational invariance and the antisymmetry of the commutator to obtain 



w = (dQ/dt) 



dtHK^KaiO)]), 



(15) 



where we have now used the brackets to indicate the expection in whatever ensemble we choose for the iV, $, L 
particles. We could also have obtained (| 1 51) starting from the standard "reduction formulae" of field theory. 
Choosing an equilibrium ensemble for N,Q,L we get the rate expression, 



w = Z _1 Tr 



dt [Ki(t),K 2 {0)] 



(16) 



where Z is the partition function and where we take the chemical potentials to satisfy the equilibrium condition 
fJ'N = + A*$- We specialize to the case in which there are no L's in the medium so that only the term with K\ on 
the left in the commutator in (|l(j|> enters, and we set the expectation of c^c] equal to unity, giving, 



w = g' 2 / dE f{E )r(E. } ). 



where 



r{E) = Z^Tr 



dte^ Et [N(t)^(tUN\0)mi 



(17) 



(18) 



We note that the RHS of l|18l) appears to depend on the index j. It is the independence of the N and <& energies 
on the respective particle momenta that makes the expression j independent, after the trace has been taken. With 
our normalizations all of the integrals over momenta in the perturbation expansion of l|18(l will be unity, and we need 
concern ourselves with the energy integrals only. 

We now use the methods of thermal field theory to evaluate the average in (|18fl . In thermal field theory, however, 
a Heisenberg picture operator is defined by 



0(t) = exp[i(H - n.i^i)t] Sc h exp[-i(H - 2J ^iMi)*] 



(19) 



and in (|18fl the time dependence of the operators N(t) and <jr(t) is given by the Hamiltonian alone. We can, however, 
use the time dependences of (|19|) in l|18[) provided that E in l|18f) is evaluated in the end at E = Ej — /ijv + l^<t> = Ej—fiL, 
rather than at uij. We thus replace fTT| by 



(20) 



g" / dEjfiEj)^-^), 



with the additional understanding that now the evaluation of r(E) in (|18f) proceeds according to the rules of equilibrium 
statistical mechanics. 

As a point of terminology, we note that the operators that enter in the definition of the function r(E), namely Nffi 
and AT T $, are the "currents" that enter the RHS of the field equations for the fields L and L* , divided by g. Thus 
the thermal average of the unordered product, as we find it in H18|) . is directly related to the imaginary part of the 
self-energy function for the L field. The construction of thermal Greens functions proceeds in imaginary time, and 
thermal factors arise in the continuation back to real energies. The result is 



r(E) = 



Im 



M(E + it]) 



(21) 
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where M(E) is the imaginary time self-energy (or "polarization part") for the L, and r\ is infinitesimal. 3 The 
argument of Mq reflects the continuation described above. Note that the imaginary part of this polarization part is 
not related to any real decay of the L, which is stable; it arises from energy-conserving processes in which the L meets 
a $ from the medium, making whatever state, before coming back to the original configuration. 

In the first part of the calculation of M(E), given below, E is assigned a value that is equal to irim/{3 where m 
is an odd positive or negative integer, according to the rule for the discrete energies of Fermions. At the end of this 
calculation we have a function of E which despite having been defined on a discrete set of points is the unique extension 
to the whole E plane that vanishes at infinite E. This function is to be continued to a nearly real energy, as indicated 
in (|2f J) , with the usual infinitesimal imaginary part in E. The quantities e$ = E$ — /!$ and ejy = En — hn remain 
real throughout this calculation. The lowest order single-loop graph for L — » L has a (N, $) loop with circulating 
energy, uj n which we take to be the energy of the $ line. The inverse propagator for the $ line is oj n — e$; then the 
energy associated with the N line is E + u> n with inverse propagator E + u> n — €n ■ The one- loop form for M, which 
we denote by Mq and which is of zero'th order in g 2 , is given by the sum over Matsubara frequencies as, 



M (^)=/3- 1 £ t wT- r, (22) 

where uj n — 2nin/0, in view of our choosing w n to be the energy of the intermediate <& (i.e. boson) line in the 
self-energy bubble. This sum can be evaluated using, 

0- 1 "If F(2*ni) = / ^icot(^V H, (23) 

n=-oo C 

which holds for a function F(z) that has no singularities on the imaginary axis. The contour, C, runs from — oo to oo 
just below the real axis, then returns to — oo just above the real axis. We get, 



For evaluation on the lower part of the contour we use, 



t)(iE- 



ie N ) 



2i V 2 



For the upper part of the contour we use, 



2i 



■ cot 



(f) 



f 

2 

1 

2 



1 



-i0u _ i ' 



(24) 



(25) 



(26) 



The contribution to Mq(E) from the ±1/2 terms in the identities 1)25(1 . (|26|l can easily be seen to vanish, since both 
poles of the multiplying factors of the integrand are in the upper-half plane. The remaining parts, coming from the 
second terms on the RHS's of (I25|l and (|26|l . are evaluated by noting that we can close the contours at infinity and 
calculate the residues at u> = ie^ and at to = — iE + iepj (remembering that E is imaginary whereas ejv and are 
real and positive), giving, 



M (E) = 



E-e N 



1 



1- 



(27) 



Note that the discrete value of the energy , E — inn/ (3 with n odd, entered here to create the Fermi distribution 
deriving from the pole at lo = —iE + itM- Next we use l(27 |l in l(2T )l . obt ammg, 



r{E) = 2T:g' 2 5(E - e N + e$)(l + e ^N-/3^yi \^ e p H _ x yi + ( e 0e N + ^-1] 



= 2Trg ,2 6{E -e N + e$)[l - e"^]- 1 ^ 



(28) 



3 The derivation of 1211 for the case of fermionic currents is given in ref llll . eq.s 4.3 to 4.8. It can also be found in ref Il2j . ch. 9, where 
it is, however, given explicitly for the Bose case, so that the appropriate change must be made in the multiplicative thermal factor in 
order to obtain 1211 . 



6 



leading finally to the rate, as given by (|20|l 

w^=2irg' 2 J dE j f{E j )8{E j +E^,-E N )[l-e-^i']- 1 [l + e^ N ]' 1 . (29) 

Eq. lj2l3|l is just the golden rule © 

for the lowest order N — » L + $ process with the appropriate thermal factors added for the N and the <I>. 

While it may seem to have been a pointless task to derive the trivial lowest order rate in the above complicated 
way, the reward is that the corrections to order g 2 are now relatively easy to define and to calculate. We need only to 
put in the correction to the N inverse propagator, replacing E -\-u> n — ejv in (|22J) by E + u> n — ejy — Tl(E + u> n ), where, 

/oo 
dErfiEj)] 1. ( 3 °) 

n — — oc ^ J 

and where £_,- = Ej — ^l- The function II is the self-energy bubble in which the N goes virtually into a $ and an L. 

There is a completely analogous correction to the $ inverse propagator, which we will not write down. If we 
made the appropriate replacement of lu h — e$ in l(22|) and expanded, we would generate the scattering process B, i.e. 
N + L^Q + L + L^N + L, plus the appropriate singularity-regulating corrections to the basic decay N — > L + $. 
The results of the calculations are completely parallel to those we will discuss that relate to process A. For simplicity 
only, we stick to the range of parameters where the B channel is disfavored. 

In the language being used in the literature about this subject area, the relevant graphs have already been "re- 
summed" when n is included in the inverse propagator. Indeed, the papers to which we compare our results keep only 
Im[II(x)], where x — ejy and the appropriate continuation has been made before taking the imaginary part. However, 
as we see below, no "resumming" of any kind is needed in order to get the complete and finite answer to order g 2 . 
Thus we simply use only the order g 2 term in the expansion of the propagator, i.e. we make the replacement in 1)22(1 . 

(E + uj- en)- 1 -> (E + uj - ew)" 1 + (E + uj - e N )- 2 U(E + uj). (31) 
The last term supplies the order g 2 correction to the rate, given simply by replacing Mo in H21|) by 

1 



M 1 (E)=g 2 0- 1 E f— 



uj„ - e$ {E + uj„ - e N ) 2 

n' — oc 1 

dEiftJSj) Y, "FT w\ (32) 

' u n > — e$ E + u n — u) n > - Ej J 

n — — oo J 

and then substituting in (|20|l . 

Actually, the structure of the expression H32|l probably is sufficient by itself to show that the answer is defined, 
and finite, although we will show below the explicit form that results from doing the u) n , u) n i sums, again using the 
identities (|25|l and J2EJl, followed be the continuation to (almost) real values of the energy, followed by taking the 
imaginary part. Indeed, the problem with the perturbation theoretic S-matrix approach that led to the perceived, 
and spurious, need for the "resummation" , is precisely the double pole in the integrand, [E + uj n — €n)~ 2 . In the 
S matrix approach (in the superficial treatments that have been followed) a closely related factor appears under an 
integral in a generic form {E + u> — E^ + ir])^ 1 (E + uj — E^ — iri)~ 1 , where i] is infitesimal, the integration contour 
being squeezed between the singularities on both sides. 

In (|32|) . no such problem arises in the evaluation. One way of seeing this, and the best first step before doing the 
frequency sums, is to write (E + uj n - e N )~ 2 = (d/de N ){E + uj n - e w ) _1 - 4 Then the sums are performed using 
followed by the continuation, taking the imaginary part, and finally taking the (or E n ) derivative. Doing the 
frequency sums, first for the inner sum in 1)32(1 . and then for the outer sum, is a straightforward mechanical procedure, 



4 This trivial trick is unavailable in the standard S matrix approach if we have resolved into exclusive final states since in that case the 
denominator function is not a square, but rather an absolute square. 
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and we omit the details. The answer, after we substitute in l|20|) is remarkably simple, 



1 



1 



1 — e -/9(E*-M*) e /3(B«-M«) — 1 



dE 



N 



1 



' Ei + Ea 



E 



-][f(E s ) - f(E N - E+) 



N 



3 l + e P( E 3-» L ) 
I _|_ gPiEj+E^-HN) 

I _|_ e /3(-Ejv-/Jjv) 



(33) 



This is the complete second order correction to the rate for L appearance in the model described by the original 
Hamiltonian. We note that the term with f(Ej) in the last factor is exactly the (divergent) rate calculated from the 
L + <f> — » L + $ Fcynman graph, without "resumming" . 

Before interpreting the other term, which removes the divergence, we take the limit of Boltzmann statistics (all /x's 
large and negative) 



w { b ] = fl yVG*+i»i.-E») j dE . f^.y-PEt 



L Ej + E^—Eiv 



][f{E 3 )-f{E N -E^)e 



(34) 



The integrands (after taking the continuum limit) in both 134|) and in l|33|> are perfectly finite; they do not require 
regulation, or even principal value prescriptions. Moreover, there is no assurance that is positive; the counter 
term with /(En — E<&) can easily outweigh the first term in the last factor in (|34(l . and does so for some simple choices 
of f(E). 



ALTERNATIVE APPROACH 



There are two other ways of obtaining the same results. In a careful real-time treatment as well, the inclusive rate 
w can be expressed as the imaginary part of a two-loop amplitude, where the propagators are standard real time 
finite temperature propagators. In a correct formulation (with irfs defining the integrals) all of the singularities are 
on the same side of the integration contours, so that this amplitude is non-singular. The real time approach may be 
less authoritative even if is a quicker way to get to the result. See ref. for a demonstration of its delicacy in 
some circumstances similar to the ones of this paper. A real time approach based on Keldysh (2x2 matrix) Greens 
functions should be problem-free however. 5 

The third approach goes back to the vacuum S matrix calculation based on the Feynman rules but treats in addition 
the medium-dependent wave- function renormalization (MDWFR) that is required in such an approach [15| . We begin 
from the lowest order result for the Boltzmann statistics case, but modified as follows, 



2ng u e 



12 J3{fi N -E N -A) 



E a 



A)Z 2 



(35) 



which is just the Boltzmann limit of JSJ but with an energy shift A for the N particle and a factor Z 2 , which is a 
redefinition of the coupling constant. Now we calculate the lowest order difference between this modified rate and the 
earlier lowest order result, 



d 



dE N 



f(E N - E 9 ) + [Z 2 - l)f(E N - E q 



The N propagator, with the single bubble self energy part is 

S(E) 



E - E N +g 2 I dEj 



Ei +E&-E 



(36) 



(37) 



5 In principle, the Keldysh Greens function approach of ref. ll4l should produce our results, both in the toy model and in the genuine 
leptogenesis models. The particular contribution that we have discussed is con tain ed in the two loop graph of fig. 12c of this paper, a 
self energy part for the I fields. However in the sample evaluation given in ref. Ill of the imaginary part of this graph, the internal N 
propagator is taken as an energy-independent (and point) interaction, so that none of the problems that we have addressed are present. 
If instead one uses the propagator that creates the singularities we have discussed in this paper, there will appear to be squar es o f delta 
functions in individual terms, which are characteristic of real time approaches. But in the matrix (Keldysh) approach of ref, ll4l these 
terms should cancel, once again, so that there is no need for "resumming" for purposes of regulation. 
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The identification of the energy shift A and Z — 1 is accomplished by calculating the position of the pole of the 
propagator and the associated residue, both to lowest order 6 , obtaining, 

A = -g 2 [dE, f( f j) , (38) 



Ej + E$ — En 



and 



z * - 1 - p + ' 2 / ^ (^W '- - 1 - -» 2 ? (^W- < 39) 

Now we see that we can replace our imaginary time calculation of the last section by adding the "Feynman graph" 
result, the first term on the RHS of l|34|) . to Sw^ from J3§}, using i|39|) and (|38fl . and defining the result as 
The difference between this result and l|34fl is given by, 

■f(E N -E<,)e-0 E » 



(i) 

mdwfr ' 



dE 



N 



Ej + E<$ — En 



e~^[-(3 + g^f(E N - E^\ e~^f(E N - 

Ej+E^-En (Ej+E^-En) 2 ) ' [ ' 

Thus we can get the identical result as Q34JI. with a lot less work in the case of the present model, using the methods 
of this section. There are some reasons to favor the imaginary time approach of the previous section : a.) Rather than 
depending on the cancellation of terms that are only defined formally, in a sense, the answer is never separated into 
individually divergent terms, b.) In the model at hand, it is the N propagator that serves to define the renormalization 
constants A and Zi- Since N is an unstable particle, these definitions could be questioned, c.) In the more difficult 
relativistic problem studied in it was found that the terms ±1/2 in i|25|) and l|26|) played the special role of picking 
out the "vacuum" part of the expression, a part that was independent of the temperature and the density of the 
medium, but that carries ultraviolet divergences that are to be addressed in the usual fashion. 



DISCUSSION 



We explain further why the toy model poses the right question: 

a. Momenta and relativity. We not only left out relativity, but we discarded the kinetic energies of two of our four 
kinds of particles. In a limit in which the N and $ are very massive, this is not such a bad mutilation, since the kinetic 
energy change for the heavy particle in any of our reactions would be small in any case. What if all the particles had 
finite mass and were relativistic? In the first place we would have a bigger manifold of non-trivial integration for the 
two loop processes. In (|32f) we had a sum over momenta for the L line only; in the realistic case there will be integrals 
over two three-dimensional loop momenta, and the constants En, E§ would be replaced by functions of the momenta. 
Moreover, there would be many further dependences of the summand on the loop energies, re-expressed as Matsubara 
frequencies, coming from the relativistic propagators. However, as regards to the double-pole singularity in energy 
which enters, when the innocent factors multiplying this pole are evaluated at the pole position, the resulting structure 
will recapitulate exactly the structure that we dealt with above. Thus we are assured that there is no "resummation" 
required to obtain a finite answer, again of order g' 2 g 2 , and therefore probably no real point in doing the detailed 
calculation. 

b. The initial state. The actual models being used in leptogenesis are of sufficient internal complexity to rather 
frustrate the non-specialist (such as the present author) in confronting some of the details in this regard. Thus we 
avoid any comments on the larger issues and stick with our model with four species. But to begin with our first 
arbitrary assumption, that the number of L's in the medium was zero, this was just for convenience, in order to have 
half as many terms to deal with. When, in the derivation of itTTjl we set Cic\ = 1, and c^cf = we could have set 



These MDWFR parameters are not even dependent on the medium, in a sense, in the present calculation. Note, however, that they 
would have been, if we had used the statistical factors in 1331 instead of taking the Boltzmann limit. 
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(cjcj) = 1 - n F (Ei - m), 

(4c i }=n F (E i - f i L ), (41) 

and calculated separately the terms for the creation and destruction of X's, and, if we liked, the back reaction on 
Taking the field L to be the true antiparticle manifestation of the field L and adding an equilibration process, 
(self-conjugate boson) <-> L + L, we could also watch the system move towards its true equilibrium with equal numbers 
of L and L particles and (1l — P-l = 0. But our point was rather to show how we can calculate the separate rates 
for the terms that increase and decrease the lepton number, and to show that these terms are individually free of 
the "resonance" pathology. Could one contrive an initial state in which the pathology was real? Perhaps one could, 
although the answer is not obvious. But, to put questions like this into context, we might reflect on the role of 
statistical mechanics in early universe problems. To the extent that we have equilibrium, everything is well defined 
and it is easy to retrieve the results. We can also address situations that are out of equilibrium by a little, or which 
are in thermal but not chemical equilibrium. The tools that we have described are sufficient for these cases. To go 
farther in taking systems out of equilibrium is tricky. We must hope that physics does not force us to follow the 
complete history of each degree of freedom. 

We note that the authors of Ref . @ agree with the qualitative conclusion (no resonant enhancement) that we express 
above. However, they achieve this result by subtracting twice as much from the "resonant" term than is implied by 
the replacement (JJJJ. Besides being arbitrary (in our opinion), this procedure leads to completely incorrect residual 
terms of order g 2 g' 2 as one can see by comparison with our results. If one wants to follow the graph-based (exclusive) 
approach, there appears to be no alternative to correctly treating the wave function renormalization processes as in 
the last section. 

To calculate the analogue of our g 2 g' 2 terms, but in a full relativistic model, one could follow the steps more or less 
outlined above, first reducing the question to that of the equilibrium average of a product of currents, then using the 
standard methods of statistical mechanics. We would argue, however, that one never needs these terms. If g is small, 
then they don't matter compared to the straight g' 2 term. If g is large, perturbation theory is inapplicable, and the 
answer is anybody's guess. Since we are not working in a situation in which coupling constants are exactly known, 
and since we do not seek precision results, these higher order effects must be irrelevant. 

We have followed a very computational path to arrive at our conclusions. There surely are considerations that 
would obviate the need for "toy-models" to demonstrate the existence of the perturbation expansion of inclusive rates 
that derive from formal expressions like that of Eq. I|18|) . even when the expansions for exclusive rates do not exist. 
We reemphasize that that is exactly what happened in the toy calculation; if we do exclusive calculations with graphs, 
the corrections to the N — > $ + L rate and the lowest order $ + L — > $ + L, each separately contain terms that are 
infinite in perturbation theory (without the i"/ device), but which cancel in the inclusive rate. A graphical calculus of 
the inclusive rate begins with a current at one end and another current at the other end, with instructions to take the 
imaginary part to get the rate. When graph has two or more loops 7 , then there are "cutting rules" which in principle 
give the inclusive rate, but where the individual cuts give the contributions of various processes. These rules, as well 
as the general field-theory background for our calculations are given in Ref. [lif. What we have shown is that it is 
advisable not to use these "cutting rules" , but instead to do some integrals and combine some terms before taking 
the imaginary part. There is probably an provable general result that the inclusive rate is free of the pathologies 
addressed in this paper. In an analogous case of radiative processes in zero-temperature field theory, this was proved 
long ago [l^|-|l8j. Thus, for example, bremsstrahlung divergences in exclusive amplitudes do not affect one's ability 
to compare laboratory cross-sections to completely perturbative predictions. 

The reader may note that in the present model the problems arise because of an unstable field, N. What if we 
change the masses so that the decay is forbidden? Then, of course, to the order that we have considered, there are no 
issues to address. If one were rash enough to look at perturbative corrections to these rates, however, involving more 
particles in the medium as well as corrections that would be called "radiative", all of the same issues will surface. 
Following the S matrix approach (adding up graphs and weighting with thermal factors) in this case might well give 
finite correction terms. But they will be the wrong terms, unless we do (finite) medium-dependent wave-function 
renormalization correctly. The issues raised in this paper are really not those of treatment of unstable particles, but 
rather those of how one does field theory in a medium. 



7 Note that in our terminology, the basic TV decay is "one-loop" and the scattering process is "two-loop' 
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To sum up with respect to the leptogenesis models: in our opinion there will not be a resonant enhancement of an 
inclusive rate for change in lepton number. If our toy model is deemed inadequate for making this point, then the 
author making this claim should at least find an alternative to the incorrect procedures of Eqs. © and (J5J; if these 
procedures are conceptually wrong in a simple case, then they are wrong in a more complicated case. 8 

We comment briefly on another problem in which exactly the same issues have surfaced, and in which "resumming" , 
in the sense used above, also gives correction terms that are both of the wrong order in coupling parameters and of 
the wrong sign. The physical situation is the scattering of a neutrino from a nucleon in the supernova environment. 
There is a "zero'th order" quasielastic rate which is just the (neutral current) vacuum rate, given by cross-section 
times density. But if one looks at a Feynman graph for the correction terms in which the nucleon that interacts with 
the neutrino is scattered by other nucleons, the (squared) first order correction to the amplitude has a singularity 
that is non-integrable, this coming from the kinematical region in which no energy is exchanged. The authors of Ref.s 
[2(J-[2i| have addressed this problem by putting a finite imaginary part in a denominator, and relating this imaginary 
part to a (strong) nucleon-nucleon scattering rate in the medium. 

But in a consistent perturbation development the singularity in question is cancelled by the interference between the 
second order amplitude, dead forward and elastic for the spectator particle, and the zeroth order result, disconnected 
from the spectator |24|. The interpretation of this last term is exactly that of medium-dependent wave function 
renormalization. We emphasize that while doing calculations this way, beginning with Feynman graphs, can get to 
the right answer if done carefully enough and in a framework that recognizes wave-function renormalization, the ab 
initio methods of the present paper are a more convincing way of getting the same results. Interestingly, although 
the conventional name given to the source of the difficulties in the neutrino problem (a "colinear divergence") is an 
entirely different name to that given to the source in the leptogenesis problem ( "unstable field" ) , the resolution in the 
two cases is almost identical. 

This neutrino problem can also serve to crystallize one's understanding of why the effect of adding an interaction 
that creates a new channel for producing something often gives a reduction in the rate of this production. In every 
such problem there is a sum rule that says that the integrals of momentum-space matrix elements of some weak 
currents, weighted by kinematical factors, are independent of the other interactions in the medium. When we turn 
on, say, the nucleon-nucleon interaction in the above example, the "strength" as defined by the sum rule, is rearranged. 
The strength often moves toward domains that are out of the kinematical domain allowed for the process that we 
are calculating. Thus we are not surprised that even when an inelastic channel is opened for neutrino scattering, by 
including N-N interactions, more than enough strength is taken away from the zero'th order (quasielastic) channel 
to reduce the net neutrino scattering rate. In the problem discussed in the present paper the sign of the first order 
correction depends on the function f(Ek). In the actual leptogenesis calculation, done correctly, it could well be that 
the first order changes would be a reduction of the reaction rate, not the increase that one finds in the literature. 

We cannot refrain from commenting on an irony of history. In the late 1960's in the particle physics world it was 
politically incorrect to express the belief that there is an underlying Hamiltonian. So people asked, "if you don't have a 
Hamiltonian, how do you find the partition function, Tr[exp(— )]?" . The classic paper by Dashen, Ma and Bernstein 
|25j answered exactly this question. This paper showed very generally how one can find the partition function from 
S matrix elements, for the case of particles interacting via two-body interactions. Of course, our questions go a little 
farther than that of finding a partition function, and the interaction in our problem is more complex as well, since we 
have particle production and annihilation. If we ask, "how do we answer our questions using only S matrix elements?" 
(which is surely possible), the answer must be "very carefully". The irony is that in the present era we work from 
a Hamiltonian, and that it is intrinsically clearer just to calculate the quantities of interest from this Hamiltonian 
directly, rather than going through the intermediary of an S matrix. 
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We should emphasize that our remarks have no bearing (that we know of) on the treatement of models with a number of Fermion 
species, closely spaced in mass, for example in Ref. I Hi . 
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